Casimir energy between media-separated cylinders: the scalar case 
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We derive exact expressions for the Casimir scalar interaction energy between media-separated 
eccentric dielectric cylinders and for the media-separated cylinder-plane geometry using a mode- 
summation approach. Similarly to the electromagnetic Casimir-Lifshitz interaction energy between 
fluid-separated planar plates, the force between cylinders is attractive or repulsive depending on the 
relative values of the permittivities of the three intervening media. 
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I. INTRODUCTION 

The sign of van der Waals and Casimir [l| forces be- 
tween media-separated plates can be tailored by care- 
fully choosing the frequency-dependent permittivities of 
the plates and the intervening media. For the simplest 
geometry of two planar slabs separated by a fluid, the 
Lifshitz theory [0| predicts a repulsive force when the 
permittivity of the fluid is intermediate between those of 
the two slabs for a large range of frequency; otherwise, 
the force is attractive. Such repulsive fluctuating forces 
in fluids have been measured recently In order to 
go beyond the simple parallel plate geometry several ap- 
proximation methods have been employed, including the 
proximity force approximation (PFA) the pairwise 
summation approach (PWS) [B], and dilute-limit expan- 
sions ^, Jj] . Exact approaches have been developed to 
treat complex geometries, including semi-analytical ap- 
proaches based on scattering theory [sl-flTI]. and on the 
mode-summation technique combined with the Cauchy 
theorem [l^, There are also fully numerical meth- 
ods based on Green functions worldline approaches 
[Tsl . and the combination of boundary methods tradi- 
tionally used to compute eigenvalues of the Helmholtz 
equation with the Cauchy theorem in order to perform 
the sum over modes Most of these computations 

have been performed for vacuum-separated dielectric or 
metallic plates. Fluid-separated complex geometries have 
been the subject of recent works, including a numerical 
study of Casimir repulsive forces and torques between 
fluid-separated eccentric cylinders (see also [l^ for 
the analysis of the Casimir force in a configuration con- 
sisting of an object contained inside a spheroidal cavity 
filled with a dielectric medium). 

Here we derive the exact analytical expression for the 
Casimir interaction in media-separated cylindrical config- 
urations, including eccentric cylinders and the cylinder- 
plane geometry. For simplicity, we consider in this work 
the case of a quantum real scalar field satisfying usual 
boundary conditions on the interfaces. Our approach is 
a generalization of our previous work on the electromag- 



netic Casimir interaction in cylindrical geometries with 
perfect reflectors to the case of dielectric media, and 
is based on the computation of the Casimir energy as 
a sum of the zero-point eigenfrequencies of the three- 
media geometry. We should note that this approach re- 
quires the notion of rea] zero-point energies, and therefore 
the three media should have negligible absorption in the 
whole range of frequencies relevant for the Casimir inter- 
action. The case of absorbing media can be considered 
with alternative techniques [l9l . [20| . 



II. MODE-SUMMATION APPROACH: 
ECCENTRIC CYLINDERS 

We consider a massless scalar field (j) in the presence of 
two parallel, eccentric dielectric cylinders of radii a and 
h {a <b) and length L ^ a,b separated by a fluid. The 
eccentricity (i.e. the distance between the centers of the 
cylinders) will be denoted by aS, with d a dimensionless 
number. The inner and outer cylinders have permittivity 
ei{uj) and £3(0;) respectively, and £2(0;) is the permittivity 
of the fluid. In order to enclose the system in a finite 
volume, we will include a very large cylinder of radius 
R ^ a,b, concentric with the outer cylinder. We will 
assume that the scalar field satisfies Dirichlet boundary 
conditions on the larger cylinder. The limit i? — >■ 00 will 
be taken at the end of the calculation. 

Using the translational symmetry along the z direc- 
tion, the solutions of the corresponding Klein-Gordon 
equation can be written as (j){r,z,t) — tf[Y)e~^'^*^'^^''^ , 
where r = -I- yy. For the region inside the inner cylin- 
der (0 < r < a), it is convenient to use a polar coordinate 
system r = (r, 0) centered on the inner cylinder. The 
most general solution for ip{v) in this region is 

^(r,0) =^C„J„(Air)e"^ (1) 

n 

where C„ are constants to be determined. For the region 
outside the cylinder of radius h it is convenient to use a 



2 



coordinate system r = (p, </>) centered on it. The solution 
in this region {b < p < R) has the form 



(2) 



where (7„,_D„ are constants to be determined. Finahy, 
in the region between the cyhnders where the fluid is 
located, the solution is 



71 

written in the (r, 6) coordinate system, or 



(3) 



(4) 



written in the (p, 4>) coordinate system. In these equa- 
tions we have defined A? = ei{ui)uj'^ / ~ [i — 1, 2, 3). 

Imposing the continuity of the field and its derivative 
at the interface between the inner cylinder and the fluid 
(r — a) one obtains 

C„J„(Aia) = AnJn{X2a) + BnH^l\\2a), 
AiC„4(Aia) = A2 [A„j;(A2a) + S„i/;«(A2a)] .(5) 

Therefore S„ = -AnlMMa)/ H';^\\2a)]R^n''^\ where 



1 - 



A2 Jre(Aia)J^(A2a) 
Ai JTi(A2a)J^(Aia) 



1 - 



A2 ./„(Aia)H;'''(A2a) 
^1 j;(AiQ)ffi^'(A2a) 



(6) 



Before imposing the boundary conditions at p = 6, we 
consider that at p = i? a, 6 the scalar field satisfies 
Dirichlet boundary conditions 



C„ J„(A3i?) + D^Mi^PiXsR) - 0. 



(7) 



The final results do not depend on whether we impose 
Dirichlet or Neumann boundary conditions at p = i?. 
As we know from previous calculations [T2j . in order 
to evaluate the Casimir energy it is convenient to ro- 
tate to imaginary frequencies lj — > i^, which implies 
that, in the above equation, J^iXsR) — > IniX^R) and 

Hk'\X3R) KnihR), where A, = ^7~(iOeWT^^- 
Taking into account the behavior of the modified Bessel 
functions /„ and Kn for large arguments, it follows that 
the coefficients C„ must vanish when i? — >■ 00 . An equiv- 
alent procedure would be, without enclosing the system 
in a finite volume, to consider only outgoing waves in 
Eq.©, which amounts to imposing C„ = 0. Either way, 
the boundary conditions at p — 5 have the form 



A3 



D^HL'H^b) = A„J„(A2fe) + S„i/,W(A26), 
D„H'y\X3b) = AM^2b) + Br,K^^\X2b). (8) 



It is then possible to find a relation between i?„ and A„ 



as Br, 



-An[Jn{X2b)/Hi^\X2b)]R'n 



-,(2-3) 



where 



1 - 



A2 j;(>2fc)gi''(>3b) 



1 



A2 -H"i''(A3b)g;'''(A2f)) 
^3 ff;<i)(A36)//i^'(A2b) 



The coefficients associated with the solution written 
with coordinates centered at the inner cylinder can be 
related to those centered at the outer cylinder by the use 
of the addition theorem for Bessel functions 



An = AmJn-TniX2a6), 
rn 

Bfi = S„j J„_„j(A2a5). 



(10) 



Combining equations ([5]), (|H]), and ([TU)) one obtains a 
linear, homogeneous system of equations. The solution 
of this system is non-trivial only if det[M] = 0, where 



1 - 



R 



(1-2) 



gr'n^(A26)J„(A2a) 

J„r{X2b)Hl^\X2a) 



Jn-m{>y2a5). 

(11) 

The interaction energy between the inner (media 1) 
and outer (media 3) cylinders is given by the sum of 
the zero-point eigen-energies that are the solutions to 
det[M] = 0. We shall assume that the absortion in the 
media is negligible and therefore the zero-point eigen- 
frequencies will be real. Moreover, as we are enclosing 
the system in a big cylinder of radius R, the solutions 
to det[M] — form a discrete set 0Jn{kz), n = 1, 2, 3, .... 
In order to compute the Casimir interaction energy, we 
will use Cauchy's theorem which involves an analytic 
continuation of the determinant to the complex plane. 
As the matrix elements in Ea. (|lip are functions of Ai — 
^ ei(Lij)uj^ I c? — fc^, there will be branch points at Ai = 0. 
For example, if the permittivities are described using the 
plasma model ei(w) = 1 — cjpj/cj", the branch points are 

located on the real axis at cji = iy^c^fc^ + cjp. , where 

is the plasma frequency of the medium i. The presence of 
branch points is typical for geometries with translational 
invariance, even for the case of a single cylinder (see for 
instance Refs. [HI,!!!]). 

The Casimir interaction energy is 



E\z — 



hcL 
47r 



/+00 
dkzY.^ujn{k.)-u^{kz)), (12) 
-00 „ 



where uj'^{kz) are the solutions of detMoo = with Moo 
given by Eq. ([TT|) with b replaced by 600 ^ b. 

Let us denote as Wmax the maximun between wi, a;2, 
and ws, the positions of the branch points. We split the 
sum over eigen-frequencies as 



(13) 



3 



where (w^ ) are the eigen- frequencies smaller (bigger) 



than I 



The sum over uj^ can be written as 



E 



1 /" d 
— i> dzz— ln[detM] , 



(14) 



where C is the contour shown in Fig|T] The contour C 
starts at ioo goes along the positive imaginary axis to 
the origin, circumvents the three branch points, follows 
the negative imaginary axis up to —ioo, and closes with 
a large semi-circunference. 




FIG. 1: Integration contour in the complex plane. Dots indi- 
cate the different branch points, where we have set (without 
loss of generality) u-i < UJ2 < 1^3. Crosses denotes the eigen- 
frequencies. There is a branch cut in the interval [—ujz,uj-i\ on 
the real axis. 



One can show that the contribution to the integral of 
the semi-circunference vanishes. Moreover, the contribu- 
tion of the segments above and below the real axis gives 
minus the sum over u}^{kz). Therefore, we end up with 
a representation of the interaction energy as an integral 
on the imaginary axis 



£^13 — 



hcL 
8^ 



dfc. 



dzz— ln[dct M/ dot (Moo)]. 



(15) 

From this point onwards, the calculation of the ex- 
act Casimir interaction energy between media-separated 
material cylinders proceeds as in the case of vacuum- 
separated perfectly reflecting cylinders 12|. The main 

fi— 2) 

difference with that calculation is the presence of Rn 

(2 — 3) 

and Rn , both of which are equal to unity in the per- 
fectly reflecting case. Following the same procedure as in 
[l^ . the interaction energy can be written as 



where the matrix elements of A are 

^n(A2a) 



Ri 



1-2) 



K,n{\2b) 1 



^ T (\_h\ f,(2-3) 



Im{^2h) R)t 



In- 



(17) 

{X2aS)Ip-m{>^2aS). 



The functions i?„ are the analytic continuation of the 
functions R„ to imaginary frequencies 



1 - 



1 - 



and 



i?^3)=i^^3)(fc„^o = 



1 - 



A2 -fn(Aia)j;(A2a) 
Ai /„(A2a)/;(Aia) 
A2 /„(Aia)A:;(A2a) 
Ai /4(Aia)K„(A2a) 



A2 ^(MiMAaft) 
A3 I„{X2b)K^{\3b) 
A2 g„(A3b)g;(A2b) 
A3 K^{\3b)K„{\2b) 



(18) 



(19) 



In order to derive this result we enclosed the system 
into a large cylinder, alternatively one can deal with 
the unbounded configuration following the approach de- 
scribed in Ref.ill. 



A. Attraction-repulsion crossover 

In this subsection we study the condition for the 
crossover from attractive to repulsive interaction. For 
simplicity, we concentrate ourselves on the case of con- 
centric cylinders. In the particular case of S — the 
matrix Anp becomes diagonal, and the Casimir interac- 



tion energy reduces to 



E- 



13 



hcL 
X In 



E 



1 



Rn 



/n(A2a) Kn{\2h) 

i?i2-3) Kn{\2a) In{\2b) 



(20) 



?(l-2) 



and Rn as a 
For the Rin~^^ 



We now analyze the signs of R, 

function of the three permittivities e. 
coefficients, we rewrite Eq. (jl8|) using the new variables 
X = and A12 = A1/A2, obtaining 



1 - 



1 - 



1 J„(Ai2a:)/;(a:) 

A12 In{x)I!^{\l2x) 

A12 K„(a:)/,',(Ai2x) 



(21) 



Since In{x), Kn{x), and I'nix) are always positive for 
X > 0, while K'^{x) is always negative, the denominator 

in Ea. (PT|) is positive. Therefore, the sign of i?!'""' IS 
determined by the numerator TV, that can be written as 



£^13 = ^ J dk, j d^ lndet(I - A) , 



(16) 



N = 1 



/(^) 
/(Ai2a:)' 



(22) 



4 



where f{x) 



It is possible to show [23] that f{x) 



is an increasing function of x. Therefore f{x) > f{Xi2x) 
for Ai2 < 1 and viceversa when A12 > 1. This imphes 
that -RI^"^' > when ei > €2 {X12 > 1), and ^1^"^^ < 
when ei < £2 (A12 < !)■ 

Using a similar argument it is possible to analyze the 

~ ^2 3') 

sign of i?„ as given by Ea.([T^. The numerator is 

. ~ (2—3) 

always positive, and then the sign of Rn is governed 
by the denominator D which reads 



Eq. p6|) . For the scalar field, "perfect conductivity" cor- 
responds to Dirichlet boundary conditions on the inter- 
faces. These conditions can be formally achieved for large 
values of the pcrmitivities ei and £3 of the two cylinders. 

Using the asymptotic expansions for In{x), Kn{x), and 
their derivatives, it follows that in the limit Xi ^ 00 the 

" (1 — 2) 

function i?„ ' takes the form 



R 



(1-2) 



1-^ 

Ai 



/;(A2a) K{X2a) 



/„(A2a) Kn{X2a) 



(24) 



D = 1 



(23) 



<M As 



where u = A36, A23 = A2/A3, and g{u) = -u^^ 
g(u) is an increasing function, it is easy to check that 
^1^"^^ < for e2 > ea, and > for £2 < £3- 

The sign of the Casimir interaction energy £^13 
in Eg. ([20]) is determined by the sign of the ratio 



R 



(1-2) 



When this ratio is positive (negative) 



the energy is negative (positive). From the above consid- 
erations it follows that for £1 < £2 < £3 (corresponding 
to M!:~^^ > and M?^^^ < 0), or for £1 > £2 > £3 (cor- 
responding to i?„ ^ < and i?„ > 0), the Casimir 
pressure is repulsive. For all other cases the Casimir pres- 
sure is attractive. Strictly speaking, the sign of the en- 
ergy is not enough to determine the atractive/repulsive 
character of the pressure. However, we have checked us- 
ing simple numerical evaluations that the integrand of 
Eg. ((20)) is always a monotonous function of b. 

As we will see in Section IV, when the cylinders are 
eccentric {6 7^ 0), the same conditions imply a repul- 
sive or attractive force, with the concentric configuration 
being an equilibrium situation. When £1 < £2 < £3 or 
when £1 > £2 > £3 this equilibrium is stable, and it is 
unstable in all other situations. Just as in the case of 
the electromagnetic Casimir-Lifshitz interaction between 
media-separated planar slabs, it is only necessary that 
the above mentioned inequalities between the three dif- 
ferent permittivities should hold in the relevant range 
of frequencies in Eq. (j20p . This frequency range is deter- 
mined by the geometrical parameters, particularly by the 
minimum distance between the cylinders. In cases where 
the three permittivities satisfy the inequalities needed for 
repulsion in some frequency range, but violate them in 
some other frequency range, the global sign of the force 
results from a competition between the different contri- 
butions to the integrand in Eq. ipO)) (see for further 
details). 



B. Perfect conductivity limit 



Similarly, in the limit A3 — > 00 the function R 
the form 



(2-3) 



A3 



/;(A26) j^;(A26) 

IniX2b) K^iMb) 



takes 



(25) 



In this perfect-conductors limit, the matrix elements of 
A can be written as Anp ~ A^^p — A^„p, where are 
the matrix elements corresponding to taking the perfect 
conductor (PC) hmit R^}'^^ = rS^^^ = 1 (see [H), and 



Ai 

-^n(A2a) 



A2 In{X2a 



/„(A2a) K,i{X2a) 

KrniX2b) 



{X2a6)Ip-rn{X2a5) 



A3 Kn{X2a) 

Irn{X2b) 



E 



(A26) 

/;„(A2&) K'M) 



Im{X2b) K,niX2b) 

(A2a(5)/p_m(A2a(5). 



The Casimir interaction energy is therefore E13 
AE13, where 



hcL 



dk. 



di Tr 



AA 



1 - APC • 



(26) 
(27) 



In this limit, and assuming that the £i are constants 
(i.e., no dispersion), it is possible to obtain the explicit 
dependence of the interaction energy Ai?i3 on the three 
permittivities. The key point is to note that 



A£;- 



13 



hcL 
4^ 
1 



dk-. 



d^^hCX2,a,b,S) 
. Ai 



+ r-/2(A2,a,6,5) 
A3 



(28) 



for some functions fi. We change variables in the 
integral above introducing polar coordinates (77, ip) in 
the plane (y^^/c, fcz), so that A2 = 77 and Ai — 
{^i/^2 — 1) cos^ (p + I for i = 1,3. The integral in ip 
can be computed explicitly and gives 



In this subsection we study the perfect conducting 
limit of the expression for the interaction energy of 



dip 



\/(ei/e2 - l)cos2(p-f 1 



^Infii 

£2, 



(29) 
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for ei, 63 00. Inserting this result in Eg. ([28)) we obtain 

(30) 

where the functions G,; involve integrals in the radial co- 
ordinate rj. 

In previous works 6] , the Casimir energy was evaluated 
using dilute-limit expansions. Eg. ([30)) gives the depen- 
dence with the permittivities in the opposite limit. The 
result shows that the convergence to the case of perfect 
conductivity is rather slow. 



The dependence of the Casimir energy on the permit- 
tivities in the perfect conductivity limit can also be de- 
rived from Eg. p2p . The procedure is similar to the one 
used to obtain Eg. ([SO)) , so we only quote the final result 



^Ell ^ ln(.i/e.) ^.p^^^ ^ M!^G-(a, d) , (33) 



for some functions G^^. In the next Section we will pro- 
vide numerical evaluations that confirm this behavior. 



III. CYLINDER-PLANE GEOMETRY 

In the case of perfect conductors, we have shown in 
[1^ that the cylinder-plane configuration is contained 
as a particular case of the exact formula for eccen- 
tric cylinders. In this Section we obtain the matrix 
elements for the cylinder-plane configuration from the 
media-separated cylinders Eqs. ([TO)) and ([TTf . As in 
[l^ . let us consider a cylinder of radius a above an in- 
finite plane. The permittivities are ei inside the cylin- 
der, 62 between the cylinder and the plane, and £3 in 
the region below the plane. Let us denote by H the dis- 
tance between the center of the cylinder and the plane. 
The expression for the interaction energy in the media- 
separated cylinder-plane geometry can be obtained from 
the eccentric cylinders formula Eq. (|16p taking the limit 
6/a, 5 00, keeping H/a = b/a — 6 fixed. Using the 
asymptotic limit of the Bessel functions, it is possible to 

show that the coefficient Rn~^^ — ^^"'"^^ , and therefore 

A3 — A 2 

it can be taken outside the sum in Eq. ([T71) . After this, 
and making use of the uniform expansion and the addi- 
tion theorem of Bessel functions, it can be also shown 
that fil 



« Y,K^{y + i)i^-P-n{y-i)^K^+p{2i), (31) 

m 

in the limit y — > 00. Finally, the expression for the matrix 
elements in the cylinder-plane geometry is given by 



A3 + A2if„(A2a) 



n+p 



(2A2iJ). (32) 



This generalizes the result of [2J] for vacuum-separated, 
perfectly conducting cylinder-plane to the three-media 
cylinder-plane geometry. 

The sign of the Casimir interaction energy for the 

~ (1 — 2) 

cylinder-plane geometry is ruled by the signs of Rh 
and (A3 — A2)/(A3 -I- A2). The attractive- repulsive char- 
acter of the force depends on the relative values of ei, 62 
and 63, in the same way as for the eccentric cylinders. 



IV. NUMERICAL EVALUATIONS 

In this Section, we show numerical results for the 
Casimir interaction energy between eccentric media sep- 
arated cylinders and a cylinder in front of an infinite 
plane. For simplicity, we consider the dispersion-less case 
in which Ci are constants. 

The value of the Casimir interaction energy is obtained 
by the numerical evaluation of Eg. ([TB)) . through the use 
of the different definitions of the matrix elements Anp 
depending upon the geometry considered (Egs. (|18|) and 
(1321) ). We numerically compute the Casimir interaction 
energy using a FORTRAN program which defines the 
matrix elements of A, computes the corresponding eigen- 
values, and finally performs the frequency and wavevec- 
tor integrations. The parameters used by the program 
are the dimension of the M matrix, the number of ad- 
dends corresponding to each element of the M matrix, 
the integration limits kz^aax S'^d. fmax, and the desired 
precision. In the following, we show the numerical re- 
sults obtained. 
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FIG. 2: (Color online). Exact Casimir interaction energy 
difference between the eccentric and concentric configurations 
as a function of 5 for different values of a = b/a. Curves 
with circles correspond to a = 2.0, while those without circles 
correspond to a = 2.5. In all cases ei = 2, £2 = 5 and £3 = 50, 
ordered as indicated in the legend. Energy is plotted in units 
of hcL. 
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Fig. [2] shows the interaction energy difference between 
the eccentric and concentric configurations as a function 
of the dimensionless eccentricity S. There are two sets of 
parameters used, namely for a — b/a ~ 2.0 and a = 2.5. 
We have used ei = 2, 62 = 5 and €3 — 50 in all cases, 
varying the relative order among them: ei < 62 < £3 
and ei > 62 > £3 for the positive curves, while ei < £2, 
£3 < £2 and £1 > £2, £3 > £2 for the negative ones. The 
plot clearly shows the change of the sign of the energy 
depending on the relative order among the dielectric con- 
stants, and also the unstable/stable equilibrium position 
at the concentric configuration. 

In order to analyze the dependence of the sign of the 
force in the cylinder-plane geometry on the relation be- 
tween £1 , £2 , and £3 we show in Figl3] the exact Casimir 
interaction energy a function of the minimum distance 
between the cylinder and the plane d = H — a. Again, 
we plot different orderings of the dielectric constants. It 
is easy to note that for a decreasing or increasing order- 
ing of the dielectric constants, we get a repulsive force, 
while in any other case the force is attractive. 
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FIG. 4: (Color online). Exact Casimir interaction energy dif- 
ference between the dielectric cylinder-plane and the perfect 
conductors configurations. As can be seen in Ea. (l33|l . the con- 
vergence to the perfect conductors case is ruled by a In ei / 
coefficient, making it slower than the naively expected from 
a simple 1/ei decay. In this plot we have set £2 = 1. Energy 
is plotted in units of hcL. 
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FIG. 3: (Color online). Exact Casimir interaction energy for 
the cylinder-plane configuration as a function of the minimum 
distance between the cylinder and the plane. In the upper 
curves, we show the energy for the cases in which ei < £2 < es 
or ei > £2 > £3, and the force is repulsive. The opposite 
case is shown in the lower part of the plot. In these cases 
£1 > £2, £3 > £2 or £1 < £2, £3 < £2, and therefore the force is 
attractive. Energy is plotted in units of hcL. 



V. CONCLUSIONS 

We have derived the exact expression for the scalar 
Casimir interaction energy between media-separated ec- 
centric cylinders, and obtained the cylinder-plane result 
as a particular case. In analogy with the well-known elec- 
tromagnetic Casimir-Lifshitz interaction energy between 
fluid-separated planes, our results show that the force in 
the two non-planar geometries studies can be repulsive or 
attractive depending on the relative strength of the per- 
mittivities of the three intervening media over a broad 
range of frequencies. We have presented both analytical 
and numerical calculations to prove that this is indeed 
the case regardless the value of the radii of the cylinders 
and of the eccentricity. 

We have considered the case of a quantum scalar field. 
The analogous calculation for the full electromagnetic 
field is much more involved. The reason is that, unlike 
the case of perfect reflectors, TE and TM modes do not 
decouple for finite values of the permittivities, and this 
fact introduces algebraic complications in the derivation 
of the exact formula presented in Section II. We expect 
to analyze this issue in a future publication. 



VI. ACKNOWLEDGMENTS 



Finally, in Fig. |4]we numerically check Eq. (p3l) with 
£3 -> 00. In this plot we show the difference between the 
interaction energy for the scalar dielectric example and 
the perfectly conducting function of log £1 / ^/e^ 

(different values of the distance to the plane are consid- 
ered). 



We thank J. Etcheverry and F. Intravaia for use- 
ful discussions. The work of DARD was funded by 
DARPA/MTO's Casimir Effect Enhancement program 
under DOE/NNSA Contract DE-AC52-06NA25396. 
FCL and EDM were supported by UBA, CONICET, and 
ANPCyT, Argentina. PIV acknowledges support of UN- 



7 



ESCO LOREAL For Women in Science Programme. 



[1] H. B. G. Casimir, Proc. K. Ned. Akad. Wet. 51, 793 
(1948). 

[2] E. M. Lifshitz, Sov. Phys. JETP 2, 73 (1956). 

[3] A. A. Feiler, L. Bergstrom, and M. W. Rultand, Lang- 
muir 24, 2274 (2008); J. N. Munday, F. Capasso, and 
V. A. Parsegian, Nature 457, 170 (2009); and references 
therein. 

[4] B. V. Derjaguin and I. I. Abrikosova, Sov. Phys. JETP 
3, 819 (1957); B. V. Derjaguin, Sci. Am. 203, 47 (1960). 

[5] V. A. Parsegian, Van der Waals Forces (Cambridge Uni- 
versity Press, Cambridge, 2006). 

[6] K. A. Milton, P. Parashar, and J. Wegner, Phys. Rev. 
Lett. 101, 160402 (2008); K. A. Milton, P. Parashar, 
and J. Wagner, arXiv:0811.0128 

[7] R. Golestanian, Phys. Rev. A 80, 012519 (2009). 

[8] T. Emig, N. Graham, R. L. Jaffe, and M. Kardar, Phys. 
Rev. Lett. 99, 170403 (2007); S. J. Rahi, T. Emig, N. 
Graham, R. L. JafFe, and M. Kardar, Phys. Rev. D 80, 
085021 (2009). 

[9] P. A. Maia Neto, A. Lambrecht, and S. Reynaud, Phys. 
Rev. A 78, 012115 (2008). 

[10] O. Kenneth and L Klich, Phys. Rev. B 78, 014103 (2008). 

[11] R. Messina, D. A. R. Dalvit, P. A. Maia Neto, A. Lam- 
brecht, and S. Reynaud, Phys. Rev. A 80, 022119 (2009). 

[12] D. A. R. Dalvit, F. C. Lombardo, F. D. MazziteUi, and 
R. Onofrio, Phys. Rev. A 74, 020101(R) (2006); F. D. 
MazziteUi, D. A. R. Dalvit, and F. C. Lombardo, New 
Journal of Physics 8, 240 (2006). 

[13] E.G. Lombardo, F.D. MazziteUi, and PL Villar, Phys. 
Rev. D 78, 085009 (2008). 



[14] A. W. Rodriguez, A. P. McGauIey, J. D. Joannopoulos, 
and S. G. Johnson, Phys. Rev. A 80, 012115 (2009); M. 
T. H. Reid, A. W. Rodriguez, J. White, and S. G. John- 
son, Phys. Rev. Lett. 103, 040401 (2009). 

[15] H. Gies, K. Langfeld, and L. Moyaerts, J. High Energy 
Phys. 06, 018 (2003). 

[16] E.G. Lombardo, F.D. MazziteUi, M. Vazquez, and P.L 
Villar, Phys. Rev. D 80, 065018 (2009); E.G. Lom- 
bardo, F.D. MazziteUi, and P.L Villar, Proceedings of 
QEEXT09; eds. Kimball A. Milton and Michael Bordag, 
World Scientific, Singapore (2010); arXiv:1003.2011 

[17] A. W. Rodriguez, J. N. Munday, J. D. Joarmopoulos, F. 
Capasso, D. A. R. Dalvit, and S. G. Johnson, Phys. Rev. 
Lett. 101, 190404 (2008). 

[18] S.J. Rahi and S. Zaheer, Phys. Rev. Lett. 104, 070405 
(2010). 

[19] E. S. S. Rosa, D. A. R. Dalvit, and P. W. Milonni, Phys. 

Rev. A 81, 033812 (2010); arXiv:0912.0279 
[20] E. Intravaia and C. Henkel, J. Phys. A: Math. Gen. 41, 

164018 (2008). 

[21] V.V. Nesterenko, J.Phys. A: Math. Theor. 41, 164005 
(2008). 

[22] V. V. Nesterenko and I. G. Pirozhenko, arXiv: 1003.08861 
[hep-th]. 

[23] See Eq. (3.7) in G.M. Joshi and S.K. Bissu, J. Austral. 

Math. Soc. (Series A) 50, 333 (1991). 
[24] T. Emig, R.L. Jaffe, M. Kardar, and A. Scardicchio, 

Phys. Rev. Lett. 96, 080403 (2006); M. Bordag, Phys. 

Rev. D 73, 125018 (2006). 



